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ABSTRACT 

Radio core dominance, the rest-frame ratio of core to lobe luminosity, has been widely 
used as a measure of Doppler boosting of a quasar’s radio jets and hence of the 
inclination of the central engine’s spin axis to the line of sight. However, the use 
of the radio lobe luminosity in the denominator (essentially to try and factor out 
the intrinsic power of the central engine) has been criticized and other proxies for 
the intrinsic engine power have been proposed. These include the optical continuum 
luminosity, and the luminosity of the narrow-line region. Each is plausible, but so far 
none has been shown to be clearly better than the others. In this paper we evaluate four 
different measures of core dominance using a new sample of 126 radio loud quasars, 
carefully selected to be as free as possible of orientation bias, together with high 
quality VLA images and optical spectra from the SDSS. We find that normalizing 
the radio core luminosity by the optical continuum luminosity yields a demonstrably 
superior orientation indicator. In addition, by comparing the equivalent widths of 
broad emission lines in our orientation-unbiased sample to those of sources in the 
MOJAVE program, we show that the beamed optical synchrotron emission from the 
jets is not a significant component of the optical continuum for the sources in our 
sample. We also discuss future applications of these results. 
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1 INTRODUCTION 

The notion of core dominance as an indicator of the orienta¬ 
tion of a radio source to the line of sight was first suggested 
in an important paper by Orr and Browne (1982). This pa¬ 
per was also the first to identify core-dominated quasars 
as ordinary FR II quasars viewed at a small angle to the 
parsec-scale jets. 

Orr and Browne defined the radio core dominance, R, 
as the ratio of the radio flux density of the quasar core ob¬ 
served at arcsecond resolution to the flux density of the ex¬ 
tended lobe emission. Both flux densities are appropriately 
K-corrected for frequency and redshift to reflect rest-frame 
values. It is widely accepted that the numerator of R is pro¬ 
portional to the Doppler beaming factor of the pc-scale jets, 
where the 5 is the Doppler factor and a is the spec¬ 
tral index (dehned as S„ oc The Doppler factor is 

5 = ( 7(1 — l3cos6))~^ where the jet has speed /3c, Lorentz 
factor 7 and makes an angle 6 to the line of sight. We ig¬ 
nore radiation from the backward pointing jet. The beaming 
factor is a strong function of the angle 6, so the core flux 
density indeed reflects the orientation of the jets to the line 
of sight. However, the core flux density must also depend 
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on the power of the central engine, and the unbeamed core 
luminosity function is broad (Cara & Lister 2008). The de¬ 
nominator in Orr and Browne’s expression is an attempt to 
factor out this dependence. 

While the flux density of the extended emission un¬ 
doubtedly depends on the time averaged power of the cen¬ 
tral engine, it depends on other factors as well, such as the 
gaseous environment and the history of the source; it may 
not be a good measure of the present engine power reflected 
by the core flux density. These points were made by Wills 
and Brotherton (1995), who suggested that the continuum 
optical flux density (they used K-corrected V-band magni¬ 
tudes) might serve as a better measure of the intrinsic engine 
power. They denoted the ratio of the radio core flux density 
to the optical continuum flux density (or equivalently, the 
ratio of the radio core luminosity to the optical continuum 
luminosity) as Rv- 

Wills & Brotherton pointed to supporting evidence for 
this connection in the results of Yee & Oke (1978) and 
Shuder (1981), who found a proportionality between the op¬ 
tical continuum luminosity and emission-line luminosity, and 
Rawlings & Saunders (1981), who found a proportionality 
between the jet power and the luminosity of the narrow-line 
region. 

They also showed that expected correlations with ori- 


© 2015 RAS 


2 K. Van Gorkom et al. 


entation dependent properties tended to be stronger with 
Rv than with R (see their Table 1). But the use of Rv has 
certainly not gained universal acceptance as the superior ori¬ 
entation indicator. Since 1996, citations to Orr and Browne 
(1982) still vastly outnumber those to Wills and Brotherton 
(1995). More recently some authors have used both R and 
Rv in their correlations (e.g. Barthel et al. (2000), Richards 
et al. (2001), Aars et al. (2005), Kimball et al. (2011)). In 
none of these papers is a strong conclusion drawn in favor 
of one measure over the other as an indicator of orientation. 

Certainly, both R and Rv are qualitative indicators of 
orientation to the line of sight because both contain the 
Doppler beamed radio core luminosity in the numerator. 
Which core dominance measure is superior then depends on 
which uses the better proxy for intrinsic engine power in the 
denominator. That is what we investigate in this paper. 

The layout of the rest of the paper is as follows. In sec- 
tionj^we describe the radio sample that we will use for most 
of the analysis. Its most important attribute is that it is con¬ 
structed so as to be as free as possible of orientation bias. 
In section [3 we use observations of MOJAVE blazars to es¬ 
timate to what extent the optical continuum flux is affected 
by beamed optical radiation from the jets. In section|^we in¬ 
troduce two additional possible estimators of intrinsic engine 
power, and their corresponding core dominance measures. In 
sections!^ and 1^ we compare and evaluate the four measures 
of core dominance. In section we apply the same analysis 
to the Hough-Readhead (1989) sample of lobe-dominated 
3C quasars. In sectionwe summarize our results. 

2 THE LOBE SELECTED RADIO SAMPLE 

The radio sample analyzed in this paper was constructed by 
Gobeille (2011) and consists of two parts. The first part is 
a high redshift sample of 123 radio-loud quasars with red- 
shifts in the range 2.5 ^ ^ 5.28 that form a complete flux 

limited sample (^ 70 mjy at 1.4 GHz). The area of sky cov¬ 
ered by this sample stretches from 7 to 17.5 hours in right 
ascension and from 0° to 65° in declination. This area is 
covered by the Sloan Digital Sky Survey (Abazajian et al. 
2003) so excellent optical data are available for each quasar. 
High resolution radio images made with the Karl G. Jan¬ 
sky Very Large Arraj0 are presented by Gobeille, Wardle 
& Cheung (2014). These are a combination of images from 
archival data for 43 sources and new A-array observations 
at 1.4 and 5 GHz for the remainder. 

The second part is a low redshift (z < 2.5) comparison 
sample of quasars from the same area of sky, with the same 
flux limit, and for which high resolution data were found in 
the VLA archive. We found 131 such objects and re-imaged 
all of them (Gobeille 2011). They too have optical data from 
the SDSS. They range from some of the 3CR quasars deeply 
imaged by Bridle et al. (1994) to sources with only brief 
observations. The vast majority of the LSS consists of 4C 
quasars and some Bologna quasars that were observed with 
the VLA as soon as optical identifications were made. Most 
of the observations are to be found in Ulvestad et al. (1981), 

^ The National Radio Astronomy Observatory is a facility of the 
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Hintzen, Ulvestad & Owen (1983), Owen & Puschell (1984), 
Gower & Hutchings (1984), Barthel et al. (1988), Price et 
al. (1993), Lonsdale, Barthel & Miley (1993). 

A critical point is that in no case was a source observed 
with any prior knowledge of its structure. Also, in this pa¬ 
per we are not counting radio sources or constructing, say, 
a luminosity function, so incompleteness is not necessarily 
an issue provided that we have a representative sample of 
sources. There is no indication that this is not the case. 

Every quasar in the combined sample has a clearly vis¬ 
ible compact radio component that is coincident with the 
optical position measured by SDSS. These cores are widely 
believed to be Doppler boosted pc-scale jets (e.g. Orr & 
Browne 1982, Blandford & Konigl 1979), and their contri¬ 
bution to the total flux density of the source is therefore 
strongly dependent on their inclination to the line of sight. 
Also, many sources exhibit prominent kpc-scale jets, which 
we know are also significantly Doppler boosted because they 
are (in quasars) visible only on one side of the source and 
this is invariably the side of the source pointing towards us 
(Garrington et al. 1988). From deep VLA images, their typ¬ 
ical velocities are found to be in the range 0.7c — 0.9c where 
c is the speed of light (Wardle & Aaron 1995). The more 
extended lobe emission is not significantly beamed, with an 
upper limit on the average hotspot advance speed of 0.1c 
(Scheuer 1995). 

If we wish to construct a sample of radio sources that is 
as free as possible of orientation bias, then we must ignore 
the emission from the cores and the jets, and include only 
those sources whose extended emission (mainly the lobes) 
by itself exceeds 70 mJy at 1.4 GHz. This was determined 
as follows. 

In order to separate the core, jet and extended emission 
flux densities, each source was treated individually depend¬ 
ing on its angular size and on the data sets available in the 
NRAO archives. 

Over 90% of the sources in the LSS have archival ob¬ 
servations available in several arrays and at several frequen¬ 
cies. Where there are both A and B array observations at 
the same frequency, they have been combined, following the 
procedure described in Murphy, Browne and Perley (1993) 
that allows for variation of the core flux density between 
observations. 

In general, the lobe flux density was calculated as the to¬ 
tal flux density minus the core flux density minus the jet flrrx 
density. These flux densities were determined in DIFMAP 
(Shepherd 1997). For the total flux density we placed a box 
around the whole source and added up all the emission. To 
guard against missing flrrx due to resolution effects, this 
was done using lower resolution observations where avail¬ 
able. The total fluxes measured in the image plane were 
also checked for consistency with the fringe amplitudes on 
the shortest baselines. 

The cores and jets were also measured in the image 
plane, by model fitting the features in DIFMAP. Often it 
was necessary to use higher frequency images to separate 
clearly the jets and cores from the extended emission. In 
particular, for sources smaller than about 10 arcsec, it was 
better to use C band or even X band images. To refer fluxes 
back to 1.4 GHz, we assumed spectral indices of 0.0 for the 
cores and 0.6 for the jets (Bridle et al 1994). 

From our combined high and low redshift samples, 126 
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Table 1. Ratio of LSS to MM broad line equivalent widths, with 
standard deviations calculated from those of the sample means 


Emission Line 

(EW>lss/(EW>mm 

H/3 

1.9 ±0.3 

Mgll 

1.6 ±0.2 

CIV 

2.7 ± 0.4 

Weighted Mean 

1.9 ±0.2 


quasars had lobe emission that by itelf exceeded 70 mjy at 
1.4 GHz, and they form what we call the “lobe selected sam¬ 
ple” or LSS. This sample is similar to the “lobe dominated 
sample” of SCR quasars defined by Hough and Readhead 
(1989) and we have six members in common. But our sam¬ 
ple is much larger and extends to much higher redshifts. It 
should be distributed uniformly in cos 0 up to a cut-off angle 
of 45° that is believed to separate FR H quasars from FR H 
radio galaxies (Barthel 1989; Wardle & Aaron 1997). 

The radio data measured from our images and used in 
this paper consist of radio core flux density at 5 GHz, the 
lobe flux density at 5 GHz (again, computed as total flux 
minus core flux minus jet flux), the projected arm length 
(defined as the angular distance from the core to the brighter 
hotspot), and the bending angle (defined as the external 
angle between two lines drawn from either hotspot to the 
core). These data are listed in Rauch et al. (in preparation). 


3 BLAZARS AND BEAMING 

If the optical continuum luminosity is to function as a proxy 
for the intrinsic engine power, it cannot show signs of signif¬ 
icant beaming. To quantify the amount of optical beaming 
present in the LSS, we compare it to a sample of blazars. The 
optical flux of a blazar is polarized and variable and is there¬ 
fore known to include a beamed synchrotron component. 
Our blazar sample is drawn from the MOJAVE-Metsahovi 
(henceforth MM) sample of 62 blazars (Savolainen et al. 
2010). We exclude sources classified as BL Lac objects leav¬ 
ing a sample of 48 highly beamed quasars. All of them have 
prominent cores and the vast majority exhibit superluminal 
motion (Lister et al. 2009). For all 48 sources, long-term total 
flux density monitoring at 22 and 37 GHz at the Metsahovi 
Radio Observatory allowed Doppler factors to be calculated 
from the flux variability (Hovatta et al. 2009). Using optical 
data for this sample from Torrealba et al. (2012), we have 
rest-frame equivalent widths (EW) for the broad lines of H/3 
(11 quasars), Mgll (24 quasars), and GIV (14 quasars). In 
Figure we show the mean rest-frame equivalent widths for 
these three broad lines for the LSS and the MM sources. 
As expected, the highly beamed quasars show equivalent 
widths that are systematically smaller than those of the LSS 
sources, consistent with the presence of beaming in the op¬ 
tical continua of the MM sources. The ratios are listed in 
Table [T] 

We expect the mean MM source to be orientated at 
a small angle to the line of sight. Thus, the luminosity of 
the optical continuum will arise from some combination of 
thermal radiation from the accretion disk and beamed syn¬ 
chrotron radiation from the jet: L,.,continuum = Ln,thermal + 



Figure 1. The mean equivalent widths of the broad lines H/3, 
Mgll, and CIV, plotted with sample standard deviations for 
the LSS sample (circles) and the MM highly beamed quasars 
(squares). The MM sources show smaller EWs across all lines, 
consistent with the presence of beaming. 


Li/,beamed(0). The mean LSS source, however, will be ori¬ 
ented at a much larger angle to the line of sight and the 
beamed contribution to the optical continuum will be neg- 
ligible. Lt,^continuum ~ Lnyhermal. 

We write the equivalent width as EW = 
t-'line! En,continuum where = / En.iinedi^. The fact that 

line 

the average equivalent width for the MM sources is about 
half that for the LSS sources implies that in the blazar sam¬ 
ple En,beamed(^) ~ E,y^thermal. But En^beamed(^) OC (5 (^)) 

where 5 is the Doppler factor and a is the spectral index. 
The median Doppler factor for the MM sources is 15.5. If 
the median source is oriented near the angle that maximizes 
superluminal motion, given by cosO = P, then 0 « 3?7. A 
straightforward calculation finds that at angles greater than 
7?1, beamed drops below 0.1L„,thermal and can be ignored. 
The fraction of sources oriented between 0° and 7?1 in 
a sample that is unbiased in orientation (i.e. distributed 
uniformly in cos 9) up to the 45° cut-off, is 2.6%. In our 
LSS sample of 126 sources, we expect about three sources 
to have 7%,beamed > 0.ILj^^thermai■ These would presumably 
be among the more core-dominated sources in the sample, 
but clearly beamed optical emission is not a major concern 
for our lobe-selected sample. 

In order to directly compare the mean equivalent widths 
of the LSS and MM samples, we should make sure that 
the results are not biased by the Baldwin effect (Baldwin 
1977), in which the equivalent widths of broad emission-lines 
vary roughly as E„ continuum- We have checked the values of 
L " for the sources in both samples and find that they 

i/,2500A “ 

cover very similar ranges. We conclude that the Baldwin ef¬ 
fect (which is quite apparent in Mgll, for which we have the 
most data) should affect both samples to a similar extent, 
and that the differences between the samples evident in Ta¬ 
ble and in Figure can indeed be attributed to beaming. 

It is, of course, not a new result that part of the optical 
continuum is beamed in some classes of quasar. For instance. 
Wills & Browne (1986) found that the [OHI] and H/3 equiva¬ 
lent widths were inversely correlated with R, the radio core 
dominance. Browne & Murphy (1987) found that quasars 
with flat radio spectra have smaller emission-line equivalent 
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widths than do steep-spectrum quasars. Wills et al. (1992) 
found a highly significant correlation between the fractional 
optical polarization and the radio core dominance, R. All of 
these results are consistent with the simple picture discussed 
above. What we have done here is attempt to quantify the 
effect of beaming, and show that it is not important in a 
sample that has been selected to be largely free of orienta¬ 
tion bias. 


4 FOUR MEASURES OF CORE DOMINANCE 

So far we have considered two definitions of core dominance, 
or equivalently, two possible proxies for the intrinsic power 
of the central engine. We can also use the Rawlings & Saun¬ 
ders (1981) result that there is a proportionality between 
the jet power and the luminosity of the narrow-line region; 
i.e. we can construct a third possible core dominance mea¬ 
sure that normalizes the core luminosity directly by the total 
luminosity of the narrow-line region. Following the prescrip¬ 
tion of Rawlings & Saunders, we take this to be Lnlr = 
3(3Z/[oii] + 1.5I/[oiii]). Where only one line is available, the 
other is calculated by the relationship I/[oiii] = 4L[oii]. 

Finally, Willott et al. (1999) relate the engine power 
to the 151 MHz luminosity by the formula Qox ~ 
3 X W, where L151 is given in units of 

10^® W Hz“^ sr“^. Punsley (2005) derived a similar expres¬ 
sion based on different physical assumptions, and referred 
to the expression derived by Willott et al. as “the Oxford 
formula.” We will do likewise, hence Qox. Given the steep 
spectra of radio lobes and the flat spectra of core compo¬ 
nents, the total luminosity at 151 MHz essentially measures 
the lobe luminosity. It has the possible advantage that it 
includes diffuse emission that is over-resolved by the VLA. 
Willott et al. make various approximations in deriving their 
expression. For instance, they assign all sources the same ori¬ 
entation and the same axial ratio. However it did not seem 
profitable to pursue more detailed models because they also, 
of necessity, assumed the same gaseous environment for each 
source (about which little is known), and the same param¬ 
eters that enter the equipartition calculation (ditto). These 
uncertainties will add source to source scatter to the cal¬ 
culation but should not change signihcantly its functional 
form. 

Here we will consider four different measures of core 


dominance: 


R = 

-^'u,5GHz Core 

(1) 

-^i/,5GHz Lobe 
-^^iVjSGHz Core 

Rv = 

(2) 

L 

u,2500A 

Rnlr = 

-£'i/',5GHz Core 

(3) 

-f/NLR 

-£/i/',5GHz Core 

II 

0 

(4) 

Qox 


Luminosities were computed from the observed flux 
densities using the concordance cosmology (Spergel et al. 
2003), and K-corrected to their rest-frame values. The 5 GHz 
core and lobe flux densities were available for all 126 sources 
of the LSS, from our VLA images. The 2500 A optical con¬ 
tinuum luminosity was taken from the SDSS quasar catalog 
(Shen et al. 2011) and was available for 99 sources. (We have 


Table 2. Spearman rank correlations among power proxies 


Parameter 

-^i/,5GHz Lobe 

L 

u,2500A 

QOx 

Lnlr 

-^UjSGHz Lobe 

1 





0 




L 

0.55 

1 



u,2500A 






< 10 -® 

0 



QOyi 

0.76 

0.54 

1 



< 10-'^ 

0.0018 

0 


Lnlr 

0.39 

0.41 

0.19 

1 


< 10-® 

< 10“® 

0.46 

0 


The Spearman coefficient p is given above the two-tailed 
probability of the correlation arising by chance. 



Figure 2. The correlation between geometric rank rgeom and 
pseudo-angle Spseudo for the LDQ sample (p = 0.53, p-val = 
0.0062). 


not corrected the optical luminosities for reddening, though 
this is not an insignificant effect. The fraction of radio-loud 
quasars that are signihcantly reddened is a function of red- 
shift and faintness, and varies from a few percent to at least 
17% (Glickman et al. 2013)). The [OH] and [OHI] line lumi¬ 
nosities for the calculation of Lnlr were also taken from the 
SDSS and were available for 66 sources. The 151 MHz hux 
densities for the calculation of Qox were taken from the 6C 
and 7G surveys (Hales et al. 1990 and Waldram et al. 1996) 
and were available for 37 sources. 

The Spearman rank correlation test shows that each of 
the four core dominance measures is highly correlated with 
all the others. This is expected because they all have the 
same numerator. The quality of each as a measure of ori¬ 
entation will depend on how well the denominator acts to 
factor out the intrinsic engine power. We attempt to eluci¬ 
date this in the next sections. 


5 CORRELATIONS AND GEOMETRICAL 
RANK 

We list the correlations among the four proxies for engine 
power in Table The correlations are significant in all but 
one case: Lnlr is not significantly correlated with Qox at 
least in this data set (p = 0.19, p — value = 0.46). This 
suggests that all four proxies do depend on engine power, but 
to differing degrees, and they are affected by other factors 
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Table 3. Spearman rank correlations with geometric rank rgeom 



Spearman p 

Two-tailed p-value 

variance 

R 

-0.27 

0.0021 

0.81 

Rv 

-0.31 

0.0018 

0.51 

Rqx 

-0.33 

0.043 

0.73 

Rnlr 

-0.25 

0.045 

0.87 


too. All of them will therefore add scatter when used to 
normalize the radio core luminosity. The “best” proxy for 
the intrinsic engine power should be the one that adds the 
least scatter. 

We now turn to geometrical arguments to help choose 
the best proxy. Although we have no a prion knowledge of 
the orientation of any radio source, information from the 
radio images can significantly constrain their orientations. 
For example, sources with the largest projected arm lengths 
(projected distance from the core to the brighter hotspot) 
must be oriented near the cut-off angle (45°). A source with 
a small arm length, however, could be either intrinsically 
large and oriented at a small angle or intrinsically small and 
at a larger angle, so its orientation is not constrained by 
this argument. Projected size therefore serves as a useful 
orientation constraint for large sources, less so for medium 
size sources, and not at all for small sources. (Our results are 
not very sensitive to the precise value of the cut-off angle, so 
we do not consider here the “receding torus” model, where 
the cut-off angle depends on luminosity (e.g. Hill et al. 1995; 
Arshakian 2005).) 

Another geometric proxy for orientation is the bend¬ 
ing angle of a source. The bending angle is here defined as 
the exterior angle produced at the intersection of two lines 
drawn from either hotspot to the core. Any misalignment 
between the two jets will be exaggerated in a source that 
makes a small angle to the line of sight. Analysis of the ob¬ 
served bending angles in the LSS shows that they are con¬ 
sistent with modest intrinsic bends drawn randomly from a 
flat distribution between 0° and 20° (Rauch 2013). Thus, 
sources with large observed bending angles 20° must be 
oriented at small angles to the line of sight. A source that 
appears straight, on the other hand, may be either intrin¬ 
sically straight or bent in the plane containing the jet and 
the line of sight, and its orientation is not constrained. 

While both geometric arguments are of limited utility 
by themselves as orientation indicators, combining them into 
a single measure usefully begins to sort the sources in the 
LSS by their orientation. We therefore combine the infor¬ 
mation in each of these geometrical measures by a summed 
rank, rgeom = r;^ -|- tba- Here, Ib is the projected linear 
separation between the core and brighter hotspot. Follow¬ 
ing Rauch (2013), we remove the observed dependence of 
source size on redshift by the empirically determined expres¬ 
sion Ib = + ^)' ^6 ussigu tho rank to each 

source sorted on ascending Ib and the rank vba on descend¬ 
ing bending angle. The geometric rank rgeom = -|- tba 

should therefore tend to sort the sources by orientation an¬ 
gle, to a considerable extent at large and small ranks, and 
less so for mid-ranks. 

Aars et al. (2005) defined a similar orientation measure 
for the Hough-Readhead (1989) sample of lobe-dominated 


3CR quasars (see section]^ that they called ^pseudo- They 
constructed it from a rank ordering of linear size and core 
dominance R. Since our goal is to evaluate core dominance 
as an orientation indicator, we cannot employ their Spaeudo in 
our analysis. In Figure]^ we plot rgeom versus ^pseudo for the 
Hough-Readhead sample. The two measures are significantly 
correlated (Spearman p — 0.53, p-val = 0.0062), which lends 
credence to the use of rgeom as a proxy for orientation in the 
present work. A useful feature of rgeom is that it is based on 
purely geometrical arguments, and is therefore essentially 
independent of radio source models or theory. 


6 TESTING THE CORE DOMINANCE 

MEASURES AGAINST GEOMETRIC RANK 


In Figure]^ we plot each core dominance parameter against 
rgeom. The Spearman correlation coefficients and their p- 
values are listed in Table The correlations of R and Ry 
with rgeom are highly significant, but those of Rox and Rnlr 
with rgeom are only marginally so. In Figure we also fitted 
a regression line to the data and calculated the variance of 
the data points about that the line. The variances are also 
listed in Tabled 

We can now address a possible concern that the (un¬ 
beamed) optical continuum emission is not isotropic. For 
instance, if it is emitted by a disk, then projection would in¬ 
troduce a cos 9 dependence in the luminosity. Also, patchy 
obscuration by the edge of the dusty torus might also de¬ 
crease the continuum luminosity at larger values of 9 (see 
DiPompeo et al. (2014) for a discussion of these and other ef¬ 
fects). That they are not large factors in the present sample 
is shown by the top two panels in Figure Decreasing con¬ 
tinuum luminosity with increasing 9 would flatten the slope 
of the regression line of Ry against rgeom, but it can be seen 
that the slopes for R and Ry are essentially identical. 

Note that the variance for Ry is significantly smaller 
than that for R. Using the Fisher F test, we find that the 
difference is significant at better than the 1% level. Thus 
Ry satisfies two of our expectations for a good orientation 
indicator: it correlates strongly with our geometric proxy for 
orientation, rgeom, and normalizing by the optical continuum 
flux introduces significantly less scatter into the correlation 
than do the other three choices. 

To explore this further we conducted the following ex¬ 
periment of randomizing the denominators. We constructed 
core dominance parameters with a shuffled power measure 
in the denominator. That is, for each source 


Rrandomized — 


Core Luminosity 

Power proxy from a different source 


(5) 


For each set of shuffled core dominances generated, we 
calculated the variance about the best fit line (corresponding 
to Figure]^. We repeated this 100,000 times for the random¬ 
izations of each of the four possible measures of core dom¬ 
inance, building up a distribution of variances . We can 
now answer to what extent any given proxy for the engine 
power is an improvement over random scatter. We plot the 
histograms of the calculated distributions in Figure]^ For 
R, Rox, and Rnlr, the values listed in Table (where the 
numerators and denominators belong to the same source) lie 
within 1(7 of the mean variance of their randomized counter- 
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Figure 3. Core dominance parameters plotted against geometric : 
left of each plot. 

parts. For Ry, however, the value listed in Tablelies 3.6 (j 
below the mean. Under the assumption that this distribu¬ 
tion is Gaussian, the probability that the value in Table 
was pulled from the randomized distribution is 0.0012. 

Figure demonstrates strikingly that normalizing a 
particular source’s core radio luminosity by the lobe lumi¬ 
nosity, the narrow-line luminosity, or the value of Qox that 
pertain to that source, is no better than normalizing by val¬ 
ues belonging to any other sources in the sample. In other 
words, they all fail as useful proxies for the intrinsic engine 
power of that source. The corresponding core dominance 
measures still have some merit as orientation indicators, but 
that stems purely from the orientation dependence of the nu¬ 
merator. The denominators mostly add noise. The exception 
is Rv where it evidently matters a great deal that the radio 
core luminosity be normalized by the optical continuum lu¬ 
minosity of the same source. We therefore conclude that Ry 
is a demonstrably superior orientation indicator, confirming 
the original suggestion made by Wills & Brotherton (1995). 

Strictly, what is demonstrated by Table [^and Figure 

is that both Lv sghz Core and L ° contain a common 
’ l/,2500A 

variable that cancels when we construct Ry. We suggest 
that this common variable is the underlying intrinsic engine 
power, which on physical grounds must enter both quanti¬ 
ties. 


7 HOUGH-READHEAD SAMPLE 

We have repeated this analysis for the smalller but very 
well-observed Hough-Readhead (1989) sample of 25 lobe- 
dominated quasars drawn from the revised SCR sample 
(Laing et al. 1983). The criteria for inclusion in the sam¬ 
ple are a FR II (Fanaroff & Riley 1974) radio morphology. 




rgeom- The variance about the best fit line is given in the lower 


and a flux density of the extended emission (presumed to be 
unbeamed) > 10 Jy at 178 MHz. 

We construct the rank rgeom as before. Here, how¬ 
ever, we adopt slightly modified definitions of projected arm 
length and bending angle, following Hough & Readhead. We 
rank on projected total linear size L (adjusted for redshift as 
before), and curvature Cj,, the angle produced at the inter¬ 
section of two lines connecting the points of brightest emis¬ 
sion on either side of the source and the core. 

We limit our analysis to R and Ry. We calculate the 
scatter about the best fit line and find it to be 0.407 in the 
case of Ry and 0.442 in the case of R. Here, the correlation 
of rgeom with R is insigificant (p = —0.22, p —value = 0.29), 
while the correlation with Ry is significant (p = —0.46, p — 
value = 0.021). As before, we generate randomized dis¬ 
tributions. The variance of Ry lies 2.3cr from the mean of 
its randomized distribution, while the variance of R lies 1.6cr 
away from its own. Fitting the randomized distributions 
to a Gaussian, we find the probability that Ry is pulled 
from randomized distribution to be 0.021. For R, we find a 
probability of 0.11. 

These results are not as dramatic as for the larger LSS, 
but they are nonetheless completely consistent with it. 


8 SUMMARY 

We have considered four measures of core dominance, each 
of which normalizes the radio core luminosity by a differ¬ 
ent observed quantity that purports to serve as a proxy for 
the intrinsic power of the central engine, and tested them 
using a lobe-selected sample (LSS) of 126 quasars. We ex¬ 
amined the possibility of contamination of the optical con¬ 
tinuum flux by beaming. Finally, we discriminated between 
the competing core dominance parameters by using geomet- 
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Figure 4. Histograms of the variance about the best fit line of the randomized core dominance measures, with calculated mean fi and 
standard deviation a. The vertical lines mark the variance tr^ of the non-randomized core dominance (from Table [^. Only Ry shows a 
variance that is significantly different from its randomized distribution. 


rical information taken from the radio images. Our results 
are summarized below: 

(i) By comparing the mean equivalent widths of the 
broad-lines H/3, Mgll, and CIV for the sources of the LSS 
with those of the MOJAVE-Metsahovi blazar sample and 
using a simple model for the continuum luminosity, we show 
that the contribution of beamed synchrotron radiation from 
the relativistic jet is not in general a significant component 
of the optical continuum flux for a sample of sources that 
has been selected to be as free as possible of orientation bias. 

(ii) We define the rank rgeom to incorporate both pro¬ 
jected arm length and bending angle as geometric proxies 
for orientation angle. We find that, when plotted against 
rgeom, Rv produces significantly less scatter about the best 
fit line than do the other three definitions of core dominance. 

(iii) By a comparison to core dominance parameters with 
randomized engine power normalizations, we find that only 
Rv produces a scatter that significantly differs from its ran¬ 
domized counterpart (S.ficr). For this reason, we conclude 
that the optical continuum luminosity performs better than 
the competing proxies for engine power at normalizing the 
central engine power, rendering Ry the superior indicator 
for orientation. 

(iv) We repeat this analysis for the Hough-Readhead 
LDQ sample, limiting our analysis to R and Ry. The re¬ 
sult is consistent with the LSS result—the variance of Ry 
lies 2.3 ct from the mean of its randomized distribution, while 
that of R lies only 1.6a away. 

This is the first of a series of papers that will explore 
methods to constrain the orientation of FR II radio sources 
in general and quasars in particular. The range of applica¬ 


tions for a reliable measure of quasar orientation is large, 
as attested to by the (currently) 458 citations to Orr and 
Browne (1982). 
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